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ABSTRACT 

Many land use and land cover and hydrologic applications require the 
use of satellite imagery ^nd data. Maps ana ^regations are made from 
the data which might exist in concert with other data in a geographic 
information system. Users have basic needs for registration and 
rectification of satellite imagery related to specifying, reformatting, 
and overlaying the data. Presently, each user must accomplish these 
tasKS independetly because the present data pre-proccosing system is 
unreliable. These data are sufficient for users who must expand much 
effort in registering data. These users have requirements concerning 
projection, pixel size, resampling, and accuracy, and most would 
satisfied with data that met the standards proposed, but not consistently 
achieved, for the present system* 

APPLICATIONS 

Users of satellite imagery .Tor land use and land cover or hydrologic 
applications generally are interested in interpreting land use and land 
cover from the data to produce maps and aggregations, inputting the data 
into geographic information systems, detecting changes over time, and 
using derived data in a predictive fashion (ORI, 1979). Ancillary data 
13 iirportant in every step, and the ability to incorporate such data into 
the analysis process directly increases utility for such applications as 
in^^entory ing , managing, and planning. The hydrologic application is 
really a subset of the more general land use a J land cover applica- 
tions, Frequently, “^articular categories of land use and land cover 
affecting the hydrologic buuget (such as irrigated agricultural land) are 
mapped and utilized in hydrologic models (much like geographic informa- 
tion syste.ms) with other data sets. 

These users interested in areas at least as large as multiple counties 
make most use of the present generation of satellite imagery. We'll 
focus on one prospective user, the U.S. Geological Survey (USGS) to 
illustrate needs most other users can be expectea to have. 

USGS began a program in ^'*74 to map land use and land cover for the 
entire Nation, to digitize the resulting polygons, and to make the data 
available for use in geographic information systems. Although present 
satellite data have not been used as primary source material, data ^ om 
improved sensor packages rray prove more useful. In any case, sa»- xite 
data can be expected to play a role in identifying areas in nee^> of 
up\iating (Milazzo, 1980), and Landsat digital data have proved ^^aluahle 
in regions where other source material is lacicing (Morrissey and Ennis, 
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The Water Resources Division of USGS has a need to map specific cate- 
gories of land use and land cover as an indirect estimator of water use. 
Field measurements and ancillary data are used to determine the average 
amount of water used for each hectare of a given land use and land cover. 
The resulting data are used in combination with other layers of ata such 
as recharge and discr^arge rates and present saturated thickness of an 
aquifer to predict future saturated thickness given alternative water use 
plans. A hydrologic model is used for this application in much the same 
way that a geographic information system might be used to relate multiple 
layers of data. 


BASIC WEEDS 

Land use and land cover and hydrologic applications users have basic 
needs for registration and rectification of satellite imagery: (1) the 

users must be able to locate data when given geographic coordinates or to 
find geographic coordinates given data coordinates; (2) it must oe 
possible to reformat the data to fit given map projections; and (3) 
scenes of the same area from different times must overlay each other. 

Presentlyf each user solv^-»s these problems independently with different 
degrees of success and at varying costs. Typically the user first finds 
a set of 30 or more control points for a Landsat scene that can be 
identified on the imagery and on the largest scale maps available. 

Methods differ somewhat. Some users identify control points from a 
digital display. Others print grayscales for areas surrounding likely 
control points on a lineprinter and correlated them with topographic 
maps. Map coordinates are determined either through measur'^ment with a 
latitude/longitude ^^r Universal Tranverse Mercator (UTM) coordinates. 
Second- or third-degree polynomials are computed by analyzing the control 
points and discarding those inaccurately determined. These polyno!aials 
are used as a calibration file to reference the data. Other data can 
then be registered to the Landsat base through the calibration file even 
taough that base might not actually be itself registered to a map. 

For display of Landsat data as a map, though, additional processing is 
required, vypicaily the Landsat data will be registe'^.d to a map base by 
specifying a calibration, desired projection, si^e of final pixel, and 
resampling scheme. This precision reformatting is a computation-intensive 
procedure, and not every user is able to accomplish it. 

Utilizin'! more than one Landsat scene, eitner to make of tempora] 
data in an analysis or to detect change with time, re^^uires registration 
of one scene to another. Again there are many ways to accomplish this 
task. Some users use a modified version of the control po^nt procedure 
used to register a scene with a map. A digital displ'v or grayscales are 
used to find control points on each scene, and calibration is established 
based on a polynomial. A less labor-intensive method is to use auto- 
correlation techniques, sometimes combined with gradient identification 
techniques, on a computer. This automatically finds control points that 
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cm then be edited. In ad*lition to the labor-savina advantages, such 
systems frequently are able to find many mere control points than the 
analyst is willing to find, say 200 rather than 50. Errors are usually 
less with more control After the equations are established, one scene 
is considered the primary scene, and the other is "rubber-sheeted” to fit 
it, another intensive computation job. Though many users are interested 
in using multitemporai data for analysis, not rveryone can carry it out, 
given the procedures required today. 

PRESENT LIMITATION 

The techniques outlined to solve present registration problems have been 
around for at least 5 years. EDIPS data have been around for 2 1/2 
years. The same procedures used befote EDIPS to solve registration 
problems are as necessary today <s they were before EDIPS became avail- 
able. Th^ re was great hope that the Master Data Pri^cessor (MDP) at 
Goddard would alleviate many of the problemr^ and provide reliably 
corrected data. The Landsat data corrected with ground contiol points 
(JC?) chosen from 1 : 24 , 000-scale maps wouxd be less than one pixel for 90 
percer.;- of the pixels in a scene. The Landsat Data-Users* Handbv.j)K also 
specified that the temporal registration offsets between two Landsat 
scenes having the same path-row location would be less than 0.5 pixel for 
90 percent of the pixels in a scene (U.S. Geological Survey, 1979). 

Those expectations are n:-t being met. The V.S, Department of Agriculture 
(USDA) Compaq ?d regxs tra*' ion quality of data from the MDP with **heir 
existing techniques and found that though they could not rely on it 
presently, they could use data from the MDP as a starting point for an 
algorichm that will automat ically register their segment data to the 
Landsat data ^ Graham and Luebbe, 1981). The Digital Mapping of Irrigated 
Cropland Technique Testing project fouri it necessary to manually find 
control points for 36 scenes in 1981 'ecause results of testirfg those 
Landsat seer us from the High Plains with assessment ratings of three or 
more where an average of no closer ten pixels from their expected 

locations (Koch, 1981) . 

The same Technique Testing project found that not only was it necessary 
to find control points for scena to scene registration (multitemporai), 
but there was more rotation of one scene ^ith respect to another with 
^ata from he MDP than there had ever been before. The increase in 
rotation requires much more computer memory for the ”rubber sh^ ting" 
algoithm. On the positive side, EDIPS data seemed to ccrrelat cetter 
than the old x-formot, presumably because of the cubic convolution 
resampling. This improvement was also noted by the California Irrigated 
Lands Technique Testing project (E. Bauer, oral commun. , 1981). 

REQUIREMENTS 

The present registration of Landsat data products : ' insufficient for the 
needs i 'and use and land cover and hydrologic app*.ications users. Much 
time und tfort is expended in referencing the data by map coordinates, 
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prec?s’c*i correcting it to fit mans, and overlaying mult iter.poral scenes. 
These steps could be simplified considerably and ev?n eliminated for some 
applications it data was really available in registered form. Jser 
demands in this arer are not extreme. Ic is fair to presume that those 
users presently dissatisfied with data registered by the MDP wouJd be 
qnite happy with data chat was registered as ,.’ell as they are able tn 
obtain themselves. 

Users with the ability to perform their own geometric corrections would 
like Landsat data to come with information on image geometry (like list- 
ings of ground control points) sufficient for facilitating geometric 
corrections (OPI, 1979). This is probably .e most common demand heard. 
Users would like listings of control points. With then, they couiu 
perform their own corrections, check out the accuracies of cJ :or- 

corrected data, or use them as a base for addi/ig additional points 
warranted. 

Those users who would like system-corrected data are divided n questions 
of projection, rotation, pixel size, resampling, nnd accuracy. Let’s 
look at each question, o'^.e at a time. 

all users will not be satisfied with data in just one standard projec- 
tion, such ns the Space Obliqo^ Mercator (SOM! :: Hocine Ob-ique Mercator 
(HCM) . Users of satellite imagery for land use and land cover applica- 
tions commonly use the data in combination with maps and data 

layers. These maps can be cast on a variety of projections, such as che 
Alters Conical Equal Area, or polar Stereographic, what is needed 
is the ability to convert data from one projection to inother. USCS does 
have a package of couputer routines, in Fortran, designed to permit tne 
transformtion of coordinate p4.irs from one map projection to another 
(U.S. Geological Survey, 1931) . Data could be made available in one cr 
t /o standard projections with so£*'ware «./ail.abie to the users o convert 
CO other projections, or the liuc of available projections could be 
expanded for the standard product with the user specifying his choice. 

^*os^ ucers would pzefet land use ar.^ land cover data that are geo- 
metrically corvect°: before delivery to them to be rotated to North 
Though arguments car. made that .•■>tation is not really necessary, the 
fact remains that when working with maps, and other data levels derived 
from maps, rotated data is a pleasure and unrotated data is a oain. 

There is a dilemna wich respect to this ques* .on. When working with an 
entire Lan-Isat c^cne or multiple scenes, a rotated data set is sub- 
stantially larger than its unrotat#^d counterpart. Whei^ wot king with a 
map quadrangle, though., rotated data fit jus*" right, and substantially, 
mo *e data ?re needed to get the quadrangle c ,iage iron the 'tr:ro;,ated 
set. 3o it really depends on the appiicatior and the size of the area of 
concern. The question of rotation should really be an option for the 
user. This v<Duld satisfy all. 

Three pixel sizes seem to vie for attention. The 57-meter square 
delivered with EDIPS seems to garner no real harsh feelings. Users 
understand that 57 meters is the resampling interval along a scan line 
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and they understand the logic of choosing that for the resampled pixel* 
There are perhaps stronger arguments to be made for going to a 50-meter 
squa^re pixel since it would correspond quite nicely to UTM coordinates 
and be compati le with the cell size of the USGS land use and land cover 
data in the grid cell form* There are also those who would like to see 
6C-meter pixels which would be easier to compare with 30-meter thematic- 
mapper data* 

There has also been much disagreement over resampling* Most users who do 
their own geometric corrections tend to use simple nearest neighbor, 
most^\ because of cost considerations* Cubic convolution does appear to 
prov i e smoother looking data. It does so by increasing the variability 
of the pixel values making the data harder to compress and perhaps 
adversely affecting the accuracy of classification algorithms that 
utilize variance. Much of what resistance there is to convoluted data is 
probably due mor " to unfamiliarity with it than to other factors* 

The question of accuracy is probably the single most important issue* At 
a minimum, users would like to know what the average accuracy of a 
particular product is. This is important above all else. If it is 
reasonable to achieve only 10-pixel accuracy for a given scene, because 
many 15-minute maps had to be used, the user needs to know* He can then 
determine whether that accuracy is sufficient for his application and 
improve upon it if necessary or possible* Perhaps the greatest dis- 
appointment from the AgRISTARS evaluation of data registration was the 
poor correlation that was found between accuracy and assessment ratings 
(Graham and Luebbe, 1981) * 

There are three distinct accuracies the user must worry about. The first 
involves the referencing problem. It may be necessary to extract pixels 
that lie within polygons digitized from a map. Most users would tell you 
that for this they require root-mean-square (RMS) accuracies at least as 
good as 1 pixel. This can be achieved, however, in two stages. The 
first might reference the data to within several pixels, while the second 
mist use local fit to achieve the sub-pixel accuracy desired* 

For fitting the data to a map, the real consideration is scale. The 
National Map Accuracy standard for horizontal accuracy is 0.02 inches at 
scales smaller than 1:20,000 (Thompson, 1979)* This translates to 254 
meters at 1:500,000, 127 meters at 1:250,000, and 51 meters at 1:100,000* 
Judging from past performance it seems reasonable to expect precision- 
corrected Landsat data to fall somewhere between the standards for 
1:500,000- and 1:250, 000-scale maps. This should be quite acceptable for 
most users. 

In a geographic information system environment, the resolution of the 
coarsest unit determines the effective resolution of the entire system* 
Since land use and land cover data ^re usually at a finer resolution than 
other data layers, such as soils, problems with inaccurate registration 
of a few pixels are not especially critical on the whole. In hydrologic 
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applications over multi-state regions, land cover data from Landsat may 
be at the finest resolution of the system even though it is reported for 
areas as large as 1 square mile* Misregistration could be critical in 
particular cases, though. Identifying residential land use within a 
specified slope interval computed from digital elevation models of 
15-meter resolution would require very good, sub-pixel registration. 
Registration becomes most critical when it come to overlaying data 
collected from the same point at different times. The California 
Irrigated Lands Technology Transfer project attains a 0.2 pixel RMS error 
for multitemporal registration (Wall and others, 1981) . Error tolerance 
depends mostly on environment and repetition. Larger errors can be 
tolerated if land use and land cover features are large and if only two 
dates are being overlaid* Those same errors become intolerable when the 
features are smaller (more boundary pixels) and when more scenes are 
.being overlaid. In general, most users would be quite pleased with the 
error limits proposed for EDIPS, that is, 0.5 pixel temporal registra- 
tion offset. 

SUGGESTIONS 

Most users are puzzled by the lack of quality registration coming out of 
the MDP. This has limited some users without access to their own 
geometric correction and overlaying algorithms and has increased costs 
substantially for the others. It is hoped that a result of the Registra- 
tion and Rectificaton Workshop will be an improvement in the accuracy of 
MD? products. One critical element in any procedure to accurately 
register satellite imagery to maps is selection of ground' control points. 
Thousands of control points have been picked by Landsat analysts in the 
past. Perhaps a way can be found to accept contributions of GCP’s from 
the users into the MDP library. There they would augment those already 
in the library, increasing the number of points available for each scene 
and the reliability of registration. 

SUMMARY 

Reliable registration of satellite imagery would greatly increase the use 
of such data for land use and land cover and hydrologic applications. 
While some users can and do accomplish their own registration, most 
cannot. The inability to work with registered data results in the loss 
of the tenpoiai dimension, production of inferior map products, and 
difficulty in using ancillary data that is registered. While sampling 
strategies satisfy users from some disciplines, land use and land cover 
are commonly interested in completing "wall to wall" surveys and produc- 
ing map products. 

Land use and land co'^er and hydrologic applications users have definite 
demands for registered satellite imagery. While more variety with 
respect to cptons is desirable, products which meet the standards 
proposed for EC."'S would be welcomed. One addition that most users 
mention io a listing of control points that can be used to check accuracy 
of a particular rrodrct or as a starting point for refinement. Satisfac- 
tion of those requirements would greatly improve the utility of satellite 
imagery for xand use and land cover and hydrologic ‘applications. 
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